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• ^ Abstract 

O We present results on quasidiffusion studies in large, 3 inch diameter, 1 inch 

c/3 thick [100] high purity germanium crystals, cooled to 50 niK in the vacuum of 
a dilution refrigerator, and exposed with 59.5 keV gamma-rays from an Am- 

p! 241 calibration source. We compare data obtained in two different detector 

'""' types, with different phonon sensor area coverage, with results from a Monte 

,—1 Carlo. The Monte Carlo includes phonon quasidiffusion and the generation of 

^ phonons created by charge carriers as they are drifted across the detector by 

|. ionization readout channels. 

in PACS numbers: 72.10.Di, 85.25. Oj, 95.35. +d 

cn 

5 1 Introduction 

L| The Cryogenic Dark Matter Search llj utilizes silicon and germanium detectors to 

search for Weakly Interacting Massive Particle (WIMP) dark matter |2l [3] candi- 
dates. The silicon or germanium nuclei provide a target mass for WIMP-nucleon in- 

^ teractions. Simultaneous measurement of both phonon energy and ionization energy 

provide a powerful discriminator between electron-recoil (relatively high ionization) 
interactions and nuclear-recoil interactions (relatively low ionization). Background 
radiation primarily interacts through electron-recoils whereas a WIMP signal would 



interact through nuclear-recoils. The experiment is located in the Soudan Mine, MN, 
U.S.A. 

The most recent phase of the CDMS experiment has involved fabrication, testing 
and commissioning of large, 3 inch diameter, 1 inch thick [100] germanium crystals. 
We present results on quasidiffusion studies in which the crystal was exposed with 
59.5 keV gamma-rays from a collimated Am-241 calibration source. Prompt phonons 
are generated from electron-recoil interactions along with Luke-Neganov |1] phonons 
created by charges as they drift through the crystal via the ionization channels' 
electric field. In our Monte Carlo, phonon transport is described by quasidiffusion [3 
[6] , which includes anisotropic propagation, isotope scattering and anharmonic decay, 
until the phonons are absorbed in either the Transition Edge Sensor (TES) [7J based 
phonon channels (photolithographically patterned on the fiat crystal surfaces) or lost 
in surface interactions. The small fraction of surface covered by the TES surfaces 
and the low surface losses result in phonon pulse constants much longer than the 
TES time constants allowing for good observation of the underlying phonon physics. 

The phonon quasidiffusion measurements presented here differ from similar heat- 
pulse experiments in which the detector is partially submerged in a pumped liquid 
helium bath [8]. These crystals are situated in vacuum, which allows for reduced 
phonon losses at the surfaces. Additionally, these crystals were cooled to 50 mK in a 
dilution refrigerator. Phonon losses in the TESs and non-instrumented surfaces can 
be determined via the partitioning of energy in the phonon channels, and phonon 
propagation parameters are studied via signal timing in the phonon channels. 

2 CDMS Detectors 

The CDMS detectors are made of [001] high purity germanium crystals with [110] 
oriented fiats. The CDMS detectors are 3 inches in diameter and 1 inch thick with 
a total mass of about 607 grams and come in two styles, the mZip and a prototype 
iZip style. 

One side of the mZip detector has two concentric, photolithographically defined 
ionization channels. The outer ionization channel functions as a guard ring and is 
used to veto events located near the radial surface. The total aluminum coverage on 
this surface is ~10%, and phonons absorbed in this surface are not measured. The 
other side of the detector has four lithographically defined phonon channels as shown 
in Figure [T] and |2} Three of the channels are arranged in equal pie shaped wedges and 
the fourth channel is circular, at large radius, and also aids in determining the events 
radial location. All four of the channels have equal area coverage. The aluminum 
coverage on the phonon- channel side is 37% active and 6.1% passive. 
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Figure 1: (left) A CDMS "mZip" detector with photolithographically defined phonon 
sensors. The crystal is 3 inches in diameter and mounted in its copper housing. The 
detector is normally run with three inner channels and an outer, annular guard 
channel. Close inspection shows six pairs of wirebonds as the outer channel was 
segmented into three, individually instrumented channels. 

Figure 2: (right) A close-up view of the aluminumtungsten phonon sensors. Phonons 
interact in the pie shaped aluminum fins where they create quasiparticles that diffuse 
into the rectangular shaped TES at the center of the "pie" . The parquet pattern is 
non-instrumented aluminum which maintains a constant potential across this detec- 
tor surface. 




Figure 3: (left) A CDMS "iZip" detector with photolithographically defined pfionon 
sensors. The crystal is 3 inches in diameter and mounted in its copper housing. The 
top surface contains an outer, guard phonon sensor and two inner phonon sensors 
from which an events x position estimate can be made. The opposite face (not 
shown) has a similar channel design, but rotated 90 degrees to determine an events 
y position estimate. 

Figure 4: (right) Close-up view of the iZip phonon channel and ionization channel 
(thin lines in between the phonon sensors). The phonon channel is held at ground 
and the ionization channel is held at ~ +(— )2 V for the top (bottom) surfaces. 



The iZip detector utilizes both anode and cathode lines on the same side of the 
detector similar to a Micro-Strip Gas Chamber (MSGC) [H] as shown in Figure [3] 
and |4j Unlike an MSGC however, these is a set of anode and cathode lines on 
both sides of the detector. This ionization channel design is used to veto events 
interacting near the detector surfaces. The phonon sensor design incorporates a 
guard ring channel, two channels on the top surface, from which x position estimates 
can be made and two channels on the bottom surface, from which y position estimates 
can be made. The total iZip aluminum coverage is, compared to the mZip, reduced 
significantly to 4.8% active and 1.5% passive per side. 

Phonons are detected when they are absorbed into the aluminum fins, break 
Cooper pairs, and these Cooper pairs diffuse into the tungsten Transition Edge Sen- 
sors (TES). The TESs are superconductors, voltage biased within the superconduct- 
ing to normal transition and the introduction of excited Cooper pairs raises their 
temperature which raises their resistance. Properly instrumenting the TESs allows 
the phonon flux power to be estimated. 



3 The CDMS Detector Monte Carlo, Phonons 

The CDMS detector Monte Carlo (CDMS-DMC) models both the charge and phonon 
dynamics; these signals combined provide information about event location within 
the detector, the amount of energy deposited in the detector and the interaction 
recoil type (electron- recoil vs. nuclear recoil). Partitioning of energy between the 
inner ionization channel and the outer, guard ionization channel combined with par- 
titioning of energy between the four phonon channels provides information about 
event location. Event location provides a discriminator against surface event con- 
tamination, a class of event vetoed in low-background dark matter search operation. 
The ratio of energy in the ionization channel vs. energy in the phonon channel indi- 
cates recoil type and provides a discriminator against gamma and beta background 
radiation, additional event types vetoed in low-background operation. If and when 
a population of WIMPs is detected, the spectrum of deposited energy will provide 
information regarding the WIMP mass [10]. A proper modeling of the CDMS detec- 
tors is therefore necessary to understand and interpret the detectors signal and to 
properly veto background events. 

In the CDMS-DMC, prompt phonons are produced at the interaction point with 
2 THz energy and isotropic distribution of wave momentum and slow-transverse, fast- 
transverse and longitudinal mode density of states weighted by {vp}^'^ , where Vp is the 
phase velocity and brackets indicate an average. Their propagation is described by 
quasidiffusion, and the CDMS-DMC incorporates phonon focusing, isotope scattering 
and an isotropic approximation for anharmonic decay. The phonon phase velocity 
given by the equation 
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where p is the crystal's mass density, 
bj is the phonon frequency, 
e^ is a component of the polarization vector e, 
c^l.(,v^ is the elastic constant tensor and, 
ka is a component of the phase velocity vector k [TT] . 
Group velocity is found by solving 



^3 1 



dk 

leading to a highly anisotropic distribution of phonons in position space [12] . A finely 
spaced lookup table of phase velocities, group velocities and polarization vectors is 
generated at the beginning of the CDMS-DMC. 



Isotope scatting causes high-energy phonons to diffuse as they propagate with a 
bullc scattering rate Tj = Bv^, and individual scattering rate given by 
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Anharmonic decay dominates at early times with a scattering rate F^ = Av^ [T2] . 
The decay rate constant A is the literature is generally given averaged over all 
modes [Hj, in the CDMS-DMC we apply anharmonic decay only to the longitudinal 
mode phonons [13] (with appropriate scaling of A). After sufficient downconversion, 
the phonons can propagate balhstically to the detector surfaces. 

At the detector surfaces, they either reflect or absorbed. Phonons absorbed on 
non-instrumented walls are removed from the MC The description of absorption 
on instrumented surfaces is more detailed however. An energy of at least twice the 
aluminum superconducting gap must be deposited in the aluminum fins while the 
remainder could either remain in the film and be readout or it could be reintroduced 
into the crystal in a downconversion process. Phonons with energy below twice the 
aluminum superconducting gap cannot be readout in our phonon channels and are 
removed from the MC. 

4 The CDMS Detector Monte Carlo, Charge 

It is desirable to include a charge Monte Carlo for numerous reasons. First, the 
ionization signal, compared to the phonon signal, provides a discriminator between 
electron-recoil and nuclear-recoil events in the CDMS detectors. Electron trans- 
port is described by a mass tensor, leading to electron transport which is oblique 
to the applied field [HI [15] and is necessary to explain and interpret signals in the 
primary and guard- ring ionization channels [TB]. Second, charges drifting through 
the detector produce a population of phonons, which contribute ~25% of the total 
phonon signal. These phonons spatial, time, energy and emitted-direction distri- 
butions should therefore be properly modeled in the CDMS-DMC. Third, phonons 
created during electron-hole recombination at the surfaces contribute ~25% of the 
total phonon signal but in a low frequency, ballistic regime that is used to provide a 
surface-event discriminator. These phonons also need to be properly modeled in the 
CDMS-DMC. Germanium has an anisotropic band structure described schematically 
in Figure |5| At low field, and low temperature, germaniums energy band structure 
in such that the hole ground state is situated in the F bands [000] location and the 
electron ground state is in the L-band [111] location. Hole propagation dynamics 
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Figure 5: Germanium band structure showing the hole ground state, F band and 
electron ground state L bands. 

are relatively simple due to propagation in the F band and the isotropic energy- 
momentum dispersion relationship thoieO^) = h'^k'^/'2m. Electron propagation dy- 
namics are significantly more complicated due to the band structure and anisotropic 
energy-momentum relationship. At low fields and low temperatures, electrons are 
unable to reach sufficient energy to propagate in the F or X-bands, and are not 
considered in the DMC The electron energy-momentum dispersion relationship is 
anisotropic and given by ee;ectron(k) = ^^/2(fcj|/?^|| +kj_/m±), where the longitudinal 
and transverse mass ratio mw/mx ~19.5. 

The SCDMS detectors are generally operated with low ~1 V potential difference 
between the ionization channels and the phonon channels (which are held at ground) . 
The potential accelerates and drift the charge carriers towards the appropriate de- 
tector face. As shown in Figure [6} the holes propagate parallel to the field, in the 
+z direction. The electrons however, due to the band anisotropy, propagate in the 
four L valleys and oblique to the electric field (~33 degrees from the —z axis). 

The charge carriers cannot accelerate indefinitely via the field and eventually 
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Figure 6: Gamma-ray event in the CDMS-DMC. Holes propagate upwards in the 
r band. The electrons propagate oblique to the electric field in the four L bands. 
Charge carriers are shown as black points. Slow transverse, fast transverse and 
longitudinal phonons are shown as red, green and blue points, respectively. 
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Figure 7: (left) Phonon energy spectra of electrons and holes drifted in [001] ger- 
manium with a 1 V field pointed in the [001] direction as determined in the CDMS- 
DMC The spectra are normalized to peak heights of unity. Many of the carriers 
have frequency <0.3 THz frequency, and are ballistic. 

Figure 8: (right) Angular distribution of emitted phonon energy relative to the 
crystals [001] symmetry axis. In equilibrium, the net phonon momentum must be in 
the [001] direction, but there is a significant amount of energy back scattered into 
the direction opposite of the net drift velocity. 



scatter off of phonons to limit their velocity to the longitudinal phase velocity. Mod- 
eling of these scattering processes is complicated by anisotropy and can be handled 
either via a Herring- Vogt [T7l[16] transformation, or rejection sampling methods [T8] . 
At these low fields, the electrons (holes) emit phonons with a peak frequency of ~0.2 
(0.3) THz (shown in figures [T] and [s]) , and are ballistic. 



5 Experiments 



In the first round of experiments, a SCDMS mZip detector was exposed with an Am- 
241 gamma source to produce a population of prompt phonons and charge carriers 
within the detector. The detector ionization channels were biased at 6 V, producing 
a large population of Luke phonons. The detector was operated with the outer 
phonon channel segmented into three channels, which were individually instrumented 
and rotated 60 degrees relative to the inner three channels. In the data analysis, 
pulse starting time delays were used to determine position estimates, from which we 
selected events such that we obtained a spatially uniform distribution of calibration 
events. 
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Figure 9: Event position estimates using the partitioning of phonon energy between 
phonon sensors. Two independent estimates are made with the three inner (triangle 
shaped partition plot) or three outer (crosshair shaped partition plot) phonon chan- 
nels. (Left) event position in a CDMS mZip calibration run. (Right) CDMS-DMC 

resuhs with {AbSphonon-side, Abs^aii) = (0.1, 0.03). 

The version of the CDMS-DMC used in this analysis did not consider the oblique 
propagation of electrons and contained an analytical model for Luke phonon pro- 
duction, with longitudinal intravalley scattering. Holding all other quantities fixed, 
we varied the phonon absorption probability {Abs phonon- side) on the side of the de- 
tector with instrumented phonon channels over the values (0.01, 0.03, 0.1, 0.15, 0.2 
and 0.3). The side of the detector with instrumented ionization channels contained 
23% of the aluminum coverage as the phonon channel side and the phonon absorp- 
tion probability was scaled accordingly. We also varied the probability of a phonon 
being absorbed on the non-instrumented walls {Abs^aii) from (0.001, 0.003, 0.01, 
0.03, 0.1 and 0.3) for a total of 36 MC variations. Phonons were absorbed without 
downconversion on all surfaces. 

The partitioning of phonon energy in the different phonon channels was used 
to determine a first moment position estimate. The inner three phonon channels 
and outer three phonon channels were considered separately to produce two position 
estimates for each event. Comparisons were made between the calibration data and 
CDMS-DMC results and models with {Absphonon-side, Abs^uaii) = (0.1, 0.03), (0.15, 
0.01) and (0.15, 0.03) providing good, qualitative agreement as shown in Figure^ 
The partition plots extent matches well, but the CDMS-DMC results are not as 
diffuse as the calibration data. 

Histograms of the energy absorbed in each of the six, phonon channels were 
produced for both the calibration data and CDMS-DMC results. Better qualitative 
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Figure 10: Histograms of phonon-flux energy deposited in the six phonon chan- 
nels. (Left) CDMS niZip cahbration data. (Right) CDMS-DMC results with 

{AbSphonon-side, Ms^,au) = (0.15, 0.03). 



agreement was seen for models with {Abs 



phonon-side, Abs^^aii) = (0.15, 0.01) and (0.15, 
0.03) as shown in Figure 10 The relative location of peaks and their heights are in 
good qualitative agreement. 

Additional measurements were also made of the phonon decay time in both the 
calibration data and the CDMS-DMC results. The mZip calibration data showed 
pulse decay times of ~300 fis. Parameters that were used in the CDMS-DMC and 
bracket the calibration data decay time include absorption probabilities of 
{Absphonon-side,Abs^aii) = (0.03, 0.03), (0.1, 0.001), (0.1, 0.01) and (0.1, 0.03) with 
associated decay times of ~400, 250, 250 and 200 /xs respectively. 

The iZip style CDMS detector also underwent a calibration run and comparisons 
were made between its measured phonon decay times and compared to the CDMS- 
DMC results. The reduced aluminum coverage, compared to the mZip, resulted in 
significantly longer pulse decay times of ~900 /is. CDMS-DMC parameters that 
gave pulse decay times that bracketed the measured pulse decays included absorp- 
tion probabilities of (v46sp/,o„o„_,ide, ^^s^aw) = (0.0013, 0.01), (0.0013, 0.03), (0.0039, 
0.01), (0.013, 0.003), (0.013, 0.01) with associated decay times of ~700, 800, 700, 
900 and 700 fis. Based on the ratio of aluminum absorption coverage in the iZip 
and mZip detectors, an absorption probability of Absphonon-side,iZip corresponds to 
Absphonon-side,raZip/7-7- The iZip Absphonon-side,iZip of 0.013 is cousistcut with the 
energy partitioning and phonon decay time measurements in the mZip style detector 
indicating good agreement between calibration data and the models with only the 
aluminum coverage area tuned down. The ratio of aluminum coverage to absorption 
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probability indicates that, when incident on an aluminum fin, a phonon has a ~21% 
averaged probability of being absorbed. 

6 Conclusion 

The phonon quasidiffusion measurements presented here occurred on large 3 inch 
diameter, 1 inch thick high purity germanium crystals. These crystals were cooled 
to 50 mK in the vacuum of a dilution refrigerator. Phonon losses in the TESs and 
non-instrumented surfaces can were determined via the partitioning of energy and 
signal timing in the phonon channels. Good agreement is seen between the mZip 
and iZip CDMS-DMC models with only a corresponding adjustment in aluminum 
surface coverage. We conclude that there is a ~21% probability of phonons being 
absorbed when incident on aluminum. There was a smaller, 0.1% chance of losing 
phonons on non-instrumented surfaces. 
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